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We investigate the one-dimensional polaron problem and calculate the ground state energy and the 
effective mass. We use a real-time renormalization group method and compare our results with first 
and second order perturbation theory, with Feynman's variational principle and with the method of 
Lee, Low and Pines. 
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Introduction. The polaron has often been examined 
since Frohlich proposed the corresponding Hamiltonian 
It serves as a standard model for various problems 
involving a nonrelativistic particle moving in a scalar 
field, e.g. the interaction between nucleons and scalar 
mesons or a single electron in a solid interacting with lon- 
gitudinal optical phonons. The physical picture is that 
the particle polarizes the environment and must drag this 
polarization with it, which affects its energy and effec- 
tive mass. To describe the polaron not only in the weak 
coupling limit one has to consider it beyond perturbation 
theory. The most successful methods for determining the 
ground state energy and effective mass of the polaron are 
those of Lee, Low and Pines Pekar et al. |]||, and 
Feynman [Q. Since polaron effects have been observed 
in low-dimensional systems, the problem has also been 
studied in two dimensions p[-pO[|. 

In this paper we will examine the one-dimensional po- 
laron problem. It can be realized e.g. for a Bloch-electron 
in a one-dimensional wire or macromolecular structure. 
The excitation of an electron will be strongly influenced 
by the interaction with optical phonons |llj] . If the con- 
duction band is partially filled, one can linearize the 
electronic spectrum and the model is exactly solvable 
by using bosonization techniques [Q. However, if the 
conduction band is empty, it is necessary to consider a 
quadratic spectrum for the electron with a bare mass mo- 
This leads to the one-dimensional Frohlich Hamiltonian 
with a constant coupling to the phonons. It is our pur- 
pose to examine the ground state energy and the effec- 
tive mass of the electron for this problem. Thereby we 
will use a recently developed real-time renormalization 
group (RG) technique |l|. While this method actually 
allows nonequilibrium descriptions in the present 

paper we have only used it to determine spectral proper- 
ties. Furthermore, we will compare our results with first 
and second order perturbation theory, with Feynman's 
variational principle, and with the method of Lee, Low 
and Pines generalized to the one-dimensional case with 
finite band width. For not too large coupling constants 
we find a ground state energy near Feynman's method 
and a value for the effective mass between the result of 



Feynman and the one of Lee, Low and Pines. 

Model. The polaron problem is modelled by the 
Frohlich Hamiltonian H — Hq + Hi with 

H Q = ^ e k c{c k + ^ huj q a \ a q (!) 

k q 

Hl=Y, M k( a -q + a q) C i + q C k- ( 2 ) 
fe . q 

c\{ck) is the creation (annihilation) operator of the elec- 
tron, its energy is given by the electron dispersion 
(hk) 2 /2mo, mo being the bare mass of the electron in 
the conduction band and a q (a q ) creates (annihilates) 
a phonon with frequency uo q . Since the interaction is 
dominated by the longitudinal optical branch, we pre- 
sume dispersionless phonons, i.e. uj q = to. While the 
electron-phonon interaction coefficients M% = AI q are 
proportional to 1/q in case of the bulk polaron, the one- 
dimensional situation involves a g-independent coefficient 
M @. We define 

M=( 4moo; 

where L is the one-dimensional normalization volume. In 
analogy to the three-dimensional case a is a dimension- 
less coupling constant. In the following we choose units 
such that h = m = u> = 1. With these assumptions the 
perturbation theory produces the ground state energy 
Eg = —a and the inverse effective mass 1/m — 1 — a/2. 
RG method. In the RG we consider the S-matrix 

s _ j, e -i J dtH(t) 

— e~ i I dtH «j 1 e ^ i J dtH i(t)i ^ 

where H\(t)i is the interaction part of the Hamiltonian 
taken in the interaction picture with respect to Hq. The 
idea of the RG is to leave the S'-matrix invariant while 
successively integrating out diagrams of different time 
scales. The procedure is schematically shown in Fig. [I]. 
For a given cutoff t c in time space, we allow only for 
correlation functions of the phonons with a time scale 
t >t c . At zero temperature, the latter are given by 
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(4) with the band-width D, the error-function 



All correlation functions with time scales shorter than t c 
are accounted for by renormalized energies and coupling 
constants. A change of t c to t c + dt c is made by applying 
three steps (i) expanding the second exponential 

in (|^) and introducing normal ordering for the phonon 
operators using Wick's theorem, (ii) integrating over the 
contractions with a time scale between t c and t c + dt c 
and (iii) resumming the operators in an exponential form. 
Consequently these operators will not be limited to zero- 
phonon and one-phonon operators any more. But we 
shall see that a good approximation is achieved for not 
too large coupling constants if we neglect double or higher 
order vertex operators. One advantage of this method 
is that we only need to consider the operators of the 
electronic system, as the phonon degrees of freedom are 
integrated out in each £ c -step. For the polaron problem 
we start with t c = involving the operators as in ([!]), (j|) 
and end up with an effective Hamiltonian for t c — > oo. 
We obtain the following RG-equations 
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Here we introduced 



The (k, ^-dependence of the interaction coefficients is 
generated during the RG-flow. The second term in (||) 
is a correction term. It is due to the fact that a time 
interval connected with a contraction becomes a single 
point in time at one RG step. In the next step this leads 
to the generation of new terms which were previously 
not present. The level broadening is also included in 
(H) since all energies become complex. The terms that 
generate the double vertex operators are of fourth order 
in M. Therefore the equations (||) and (||) contain the 
order M 4 exactly. 

Feynman's method. To compare the results of the 
real-time RG we applied Feyman's method to the one- 
dimensional polaron with a finite band-width. Feynman 
used a variational principle in the path-integral formal- 
ism. Within his two-particle approximation one obtains 
for the ground state energy 
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v and w are chosen such that E g is minimum. For v = w 
the result of the perturbation theory of first order in a is 
reproduced. Following Feynman we treat small couplings 
by setting v = (1 + S)w and w = 1 |l5[| . Considering 5 
small one can now expand the right-hand-side in (Q) and 
minimize it with respect to S. The resulting energy is 

2 4 / 
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Within Feynman's approximation we obtain for the ef- 
fective mass: 
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Since (j?]) follows from a mininum principle Feynman's 
method leads to more accurate results for the energy 
than for the mass (for quantitative studies concerning 
the accuracy in three dimensions see [T^]). 

Our results for the energy using Feynman's approach 
as well as our values following from perturbation theory 
of second order are in good agreement with JlO| . 

Results. The differential equations (||) and (||) are 
solved numerically. Regarding the oscillating terms in 
these equations one recognizes that given a certain dis- 
cretization in q-space one obtains large errors with in- 
creasing frequency t c . To avoid this the phase A q k t c has 
been interpolated in q-space. 

Unfortunately (||) and (^|) do not show a convergent 
behaviour for the ground state energy for t c — > oo. One 
reason is that there are undamped modes corresponding 
to high excitations in the q-sums leading to increasing ef- 
fects on the ground state energy. In this context another 
problem arises from the fact that the correlation function 
in ((J) is not decaying. As a consequence, oscillations de- 
(7) cay as a function of t c but reoccur for sufficiently large t c 
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|L7j . The idea of our solution is to neglect further renor- 
malizaton effects of both A^ and M? for t c larger than 
a certain point tf. By doing this we obtained a damped 
oscillation of the ground state energy eo in i c -space for 
t c > t f . Therefore (§) can be integrated analytically 
which leads to 



E g = lim e (t c ) 
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In Fig. U the solution of (JT^) is shown for different 
t f . The problems mentioned above make it necessary to 
choose a finite tf where the renormalization effects be- 
yond perturbation theory are contained but the numer- 
ical instabilities do not yet occur. We choose tf = 2.5 
for all values of a. At this point only low excitations 
(A < 0.4) are not integrated out yet. Since A ~ 1 sets 
the scale for the first excited state, it is reasonable to 
assume that excited states do not have further impor- 
tant renormalization effects. The change of E g between 
tf = 2 and tf = 3 is approximately 2% for a = 0.5. 

The a-dependence of E g is shown in Fig. ||. One no- 
tices that for a < 0.7 we obtain lower values for the 
ground state energy than those of Fcynman's method. 
For larger couplings our method is no longer reliable and 
yields worse results (see Fig. ^) which is not surprising 
since we neglected double and higher order vertex oper- 
ators. 

We also calculated E g for different band- widths D. As 
one can see in Fig. |]for larger band- widths the deviation 
to the results of Feynman's method grows. 

Another quantity of interest is the mass of the polaron. 
If one differentiates (||) twice with respect to k one ob- 
tains a differential equation for 1/m. However, its solu- 
tion oscillates with increasing amplitude. Therefore for 
determining the mass we calculate both the ground state 
energy and low excited energies using the same procedure 
as above. The resulting dispersion gives us a value for the 
mass. The accuracy for the mass is worse though. For 
a = 0.3 the change of 1/m is approximately 6% between 
tf = 2 and tf = 3. 

The results for different a are shown in Fig. ^. For 
small couplings (a < 0.4), we find a value for the mass 
between the variational principle of Feynman and the 
one of Lee, Low and Pines ps| . For larger couplings 
the numerical solution is too unstable to make definite 
statements from the RG-approach. From Fig. ^ we see 
that our mass depends only slightly on the band-width. 

In summary, we used a new renormalization group 
method to study the polaron problem. Although there 
arise some difficulties from the structure of the flow- 
equations, we were able to investigate the energy be- 
yond perturbation theory. The mass was calculated as 
well, but its accuracy is worse than that of the energy. 



The results were compared to the values following from 
perturbation theory, Feynman's method and the one of 
Lee, Low and Pines. The performance of the RG is good 
for couplings a < 0.5. This restriction is due to the 
neglecting of double and higher order vertex operators 
and the undamped oscillating correlation function of the 
phonons. 
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FIG. 1. Scheme of the renormalization group. 
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FIG. 4. Ground state energy of the polaron as a function 
of D for a = 0.5. Circles: Renormalization group method de- 
scribed in this paper with tf = 2.5. Dashed and long-dashed 
line: Perturbation theory of order a and a 2 . Dotted line: 
Feynman's method. 
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FIG. 2. Ground state energy of the polaron in the renor- 
malization group method as a function oft/ for D = 10. Solid 
line: a = 0.1. Dotted line: a = 0.3. Dashed line: a = 0.5. 




RG method 

perturbation theory of first order 

perturbation theory of second order 

method of Lee, Low and Pines 

Feynman's method 



FIG. 5. Inverse mass of the polaron as a function of a for 
D = 10. Solid line: Renormalization group method described 
in this paper with tf = 2.5. Dashed and long-dashed line: 
Perturbation theory of order a and a 2 . Dotted line: Feyn- 
man's method. Dot-dashed line: result of Lee, Low and Pines. 
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FIG. 3. Ground state energy of the polaron as a function of 
a for D = 10. Solid line: Renormalization group method de- 
scribed in this paper with tf — 2.5. Dashed and long-dashed 
line: Perturbation theory of order a and a 2 . Dotted line: 
Feynman's method. 



FIG. 6. Inverse mass of the polaron as a function of D for 
a — 0.3. Circles: Renormalization group method described in 
this paper with tf = 2.5. Dashed and long-dashed line: Per- 
turbation theory of order a and a 2 . Dotted line: Feynman's 
method. Dot-dashed line: result of Lee, Low and Pines. 
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